Coherent dimer formation near Feshbach resonances in Bose-Einstein condensates 
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The results of a recent experiment with 85 Rb Bose-Einstein condensates are analyzed within the 
mean-field approximation including dissipation due to three-body recombination. The intensity 
of the dissipative term is chosen from the three-body theory for large positive scattering lengths. 
The remaining number of condensed atoms in the experiment, obtained with applied magnetic field 
pulses, were used to adjust the intensity of the dissipative term. We found that the three-body 
recombination parameter depends on the pulse rise time; i.e., for longer rise times the values found 
become consistent with the three-body theory, while for shorter pulses this coefficient is found to 
be much larger. We interpret this finding as an indication of a coherent formation of dimers. 
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Bose-Einstein condensation of dilute atomic gases has 
been observed for the first time in 87 Rb atoms [lj . Since 
then, it has been condensed atoms of 23 Na 0, 7 Li 
^ ||, metastable 4 He @, 85 Rb |, 41 K and, more 
recently, 133 Cs 8j. In Ref. |j, it was theoretically pre- 
dicted that Feshbach resonances could vary the scatter- 
ing length of atoms in systems of dilute alkali gases over 
a wide range of values. A Feshbach resonance occurs 
when the quasimolecular bound state energy is tuned 
to the energy of two colliding atoms, applying an exter- 
nal magnetic field to the system. The phenomenon was 
first realized in a Bose-Einstein condensate in Ref. [To| . 
This opened the possibility of exploring new regimes of 
BEC, allowing changes of the two-body scattering length 
from negative to positive and also from zero to infinity. 
This technique was used for the condensation and col- 
lapse control of 85 Rb atoms in the hyperfine state (F = 
2, rrifp = —2) [|| 0, 0. It was also demonstrated in 
Ref. [l3| that strongly enhanced inelastic three-body col- 
lisions occurs near Feshbach resonances. In a more recent 
experiment, in Ref. jl4j . it was explored the region of very 
large scattering lengths (up to ~ 4000 Bohr radius). 

The scattering length a has been observed to vary as 
a function of the magnetic field B, according to the the- 
oretical prediction [15|. as 

where a is the scattering length, is the background 
scattering length, B r is the resonance magnetic field of 
85 Rb, and A is the resonance width. 

In the case of 85 Rb, one has resonance width A = 11.0 
G, resonance field B r = 154.9 G and background scat- 
tering length at = — 450ao 01 1 where ao is the Bohr ra- 
dius. So, given the experimental functional dependence 
B = Bit), one can determine a = ait) and, consequently, 
the dynamics of the system in such physical conditions. 
Further, mainly in strong interaction regime, a Bose- 
Einstein condensate shows inelastic loss processes that 
cause its depletion. The dominant process of losses has 



been verified to be the three-body recombination |14lllfij . 
with a time dependence concerning a simple constant rate 
equation. All earlier observations in BEC experiments 
were consistent with a description of mean field includ- 
ing such a loss process. But, in experiments with 85 Rb 
realized in the strong interaction regime [T^j , this picture 
indicates model breakdown. Bose-Einstein condensates 
initially stable were submitted to magnetic field pulses 
carefully controled in the vicinity of 85 Rb Feshbach reso- 
nance, aiming to test the strongly interacting regime for 
diluteness parameter x — n ° 3 varying from \ — 0-01 to 
X = 0.5. The loss of atoms from BEC occurred in impres- 
sively short time scales (up to two hundreds of us) and 
disagrees with previous theoretical predictions Such 
experiments reveal higher loss of atoms in shorter mag- 
netic field pulses applied on BEC and, previously, one 
knew that as longer is the time spent near a Feshbach 
resonance as higher is the loss of atoms from BEC |(| 
(consistent with a mean field approach where the inelastic 
loss term has a constant dissipative rate). According to 
Ref. ) the results could indicate the existence of a new 
physics, that cannot be described by the Gross-Pitaevskii 
(GP) formalism. Motivated by this observed discrep- 
ancy, we investigate the dynamics of 85 Rb Bose-Einstein 
condensates submitted to such external conditions and 
time scales, when we vary the s-wave two-body scatter- 
ing length in the region of strong interaction. So, we con- 
sider a generalized mean field approach that includes the 
time dependence for both mean field coupling and three- 
body recombination parameter. Previous results consid- 
ered a mean field approach with constant mean rate for 
the three-body recombination and constant value for the 
mean field coupling [H In the present work, our 

first task is to reproduce the experimental data with time 
dependent parameters or at least verify possible limita- 
tions of the time dependent mean field approach. Next, 
we consider the magnitude of the recombination rate as 
described in literature 
with time functional dependence. 

For describing a BEC in the framework of the mean 
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field approximation, we remind that such an approach is 
valid for very dilute systems when the average inter par- 
ticle distances d are much larger than la] and the particle 
wavelengths are much larger than d jl^H^- Besides, it 
is important to pay attention to the time scales present 
in the system: the physical conditions must not change 
fast enough in order to allow the replacement of a true 
interatomic potential by the contact interaction. It is 
possible in principle for the rate of change to be larger 
than Ti/ma 2 for extremely very short changes in the in- 
teractions [2^, but it is reasonable to assume valid this 
time dependent approach at least for longer pulses. 

We begin our description from an effective Lagrangian 
of the nonconservative system as in 30], in which one 
describes the dynamics of a trapped Bose-Einstein con- 
densate in spherical symmetry with such a GPE gener- 
alization, in which one also considered losses from BEC 
by three-body recombination. This Lagrangian leads to 
the following equation of the system: 
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where u> is the geometric mean trapping frequency, Uq = 
Uo(t) — Airh 2 a(t) / m and K3 = K^{t) is the recombi- 
nation loss parameter. The wave-function '5 = ty(r,t) 
is normalized to the number of atoms N. The three- 
body recombination rate K3 is here introduced for de- 
scribing atomic losses from the condensate when three 
atoms scatter to form a molecular bound state (dimer) 
and a third atom; so, the kinetic energy of the final state 
particles allows them to escape from the trap. Other 
nonconservative processes as amplification from thermal 
cloud and dipolar relaxation are neglected, since the lat- 
ter has a much smaller effect than three-body recombi- 
nation 1 14 . Ilq and in JILA experiments |14| the thermal 
cloud is negligible (only 1000 atoms in a sample of 17500). 

A theoretical prediction of K3 is a hard task since it 
is sensitive to the detailed behavior of the interaction 
potential [25|,|2(|. However, such a calculation becomes 
simpler if we consider that a is the only important length 
scale (reasonable in the weakly bound s-wave state limit) 
and this has been considered in many works 0, 1201 121I 
HHHHI. Following Ref. , the recombination 
rate is written as 



K 3 (t) - K-[a(t)] 4 
m 



(3) 



where k should correspond to the universal value n = 
3.9. But in Ref. [H and in Ref. HJ, it was found 
< k < 65 and in Ref. < k < 67.9, that is, k 

is not universal (in general K3 depends on a three-atom 
scale|23j). In Fig. ^ we schematically give two charac- 
teristic pulses employed in JILA, in which B$, Bh, B m , 
t r and th, correspond to the initial field, the hold field, 
the final field, the rise time and the hold time of the 
pulse, respectively. For describing the dynamics of con- 
densates subjected to we consider in our calculations 
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FIG. 1: Triangular (insertion) and trapezoidal shapes of mag- 
netic field pulse (Gauss) applied to 85 Rb BEC as function of 
pulse time (/ts). The rise time t r of the pulse amounts to 
12.5 fis and the hold field is 156.7 G (scattering length a = 
2300ao). The hold time of trapezoidal pulse is 120 fis. 



the same experimental parameters and conditions used 
in |14J : as a is known to be a function of the magnetic 
field B by means of we only use time dependent 
shapes of experimental magnetic field pulses employed 
in JILA (triangular and trapezoidal pulses). The hold 
field is Bh — 156.7 G, and the end field in which one 
measures the remaining number of particles N r of the 
system is B m = 164.5 G at t = 700 /is. Further, we put 
initial field B = 166 G, corresponding to a harmonic os- 
cillator state of the system (a ~ 0), applied to an initial 
sample of iVo=16500 condensed atoms of 85 Rb. Further, 
in our approach we have used spherical symmetry with 
mean geometric frequency u> = (cjfujz) 1 ^ , for simulating 
the cylindrical geometry of JILA (radial: ui r = 2ir x 17.5 
Hz and axial: lo z = 27r x 6.8 Hz). Our time dependent 
calculations started with a Gaussian shape wave- function 
which we numerically evolve by means of Eq . (0 , using 
the Crank-Nicolson algorithm, as in Refs. I3ll l32| . We 
analyzed the loss of condensed atoms like in [l4[ , by con- 
sidering hold times from thoid = (triangular) or units 
of /is (shorter trapezoidal pulses) to nearly hundreds of 
[j,s (longer trapezoidal pulses). The behavior of the scat- 
tering length, as a function of total time of the pulse, in 
the region of strong interaction atom-atom, follows sim- 
ilarly the behavior of the employed field pulses given in 
Fig. n In triangular pulses, there is a sharp peak in the 
resonance region; and a plateau with maximum value for 
a and £ (minimum value of the field), when B = Bh is 
kept constant during the time th- 
in the upper frame of Fig. |3 we used symmetric rise 
and fall times (t r = 12.5/is, t r = 25/zs and t r = 75/zs), 
to determine the remaining fraction of atoms N r /No 
in the 85 Rb BEC as function of the hold time, by 
adjus ting our curve with the first point of Fig. [21 in 
Ref. [lj. For t r = 12.5/xs, we found k 1800, very 
far away of the values described earlier in the literature 
0mmmmmmi26|. The results show the 
same exponential decay and good concordance with ex- 
perimental data, mainly for short th (circles in upper 
frame of Fig. [21 . However, as one can realize, for longer 
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FIG. 2: Remaining fraction of atoms in 85 Rb condensate ver- 
sus hold time for some rise times (indicated inside frames). 
For upper frame, k, = 1800 and it decreases up to k ~ 100 
(lower). The initial number of particles is 16500. N r is calcu- 
lated as if it was measured at r = 700 fis (B m =164.5 G). 



hold times, experimental data point out a higher dissi- 
pation when compared with our simulations. For other 
short rise times in this frame or longer rise times (lower 
frame) , we have similar behavior, but we have to decrease 
k for better adjusting to experimental data. So, the com- 
parison with experiments show that k depends signifi- 
cantly on values of t r and th- As we know |33J, the mean 
field approach should make more adequately if we were in 
a slower process. So, we also tried to verify if our results 
would give a lower value of k (inside interval described in 
literature) if we calibrated our calculation with the last 
point of the longest pulse of JILA 0] (t r = 252.6 /is). 
Really, we found k = 100 for this case (clo ser to values 
described in literature |l9L |20L t2lL 123 . 123 . r24L I25L |2o | ^ and 
the results reproduce the experimental data for longer 
rise time but they do not make very well for shorter rise 
time, as we can observe in lower frame of Fig. [21 The 
very large value of k leads us to conclude that there is 
a coherent formation of dimers, that occurs up to nearly 
t r ~ 100/is, as one can realize observing the lower frame 
of Fig. |21 for shorter rise times. It is physically sensi- 
ble that, for shorter pulses, the presence of coherence in 
the formation of dimers would be more plausible than 
for longer pulses. So, the present calculation included 
the variation of k with the time parameters of the pulse 
and the amplification of the three-body recombination 
rate can be associated with the coherent formation of 
dimers in the inelastic collisions. Together with the sig- 
nificant loss from the coherent formation of dimers, one 



TABLE I: Numerical values of three-body recombination co- 
efficient k as function of the rise time t r of the magnetic field 
pulses applied to the 85 Rb BEC. 
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should also observe a burst of atoms carrying the excess 
of energy, which for 85 Rb with the maximum value of 
a = 4000a , would be above 70nK. Indeed, in the JILA 
experiment, it was seen a significant number loss from 
the BEC for pulses lasting only few tens of microseconds, 
which were accompanied by a burst of few thousand rel- 
atively hot (~ 150nK) atoms that remained in the trap 
|14|. In our description, for each hot atom one dimer 
is also formed. Therefore, the burst of atoms should be 
accompanied by a burst of weakly bound dimers. Using 
the observed temperature and momentum conservation 
of the recombination process, we predict that the dimers 
are also relatively hot (~ 75nK). 

We note that, if we consider a constant value of k, 
in all the cases (for any choices of th) we obtain a de- 
creasing behavior of N r , as we increase the rise time t r . 
This is in contrast with the experimental results. So, 
we consider to adjust the values of k that approximately 
better describe the experimental data for each given rise 
time t r ; i.e., k is taken as a function of the rise time 
(k = K,(t r )). Our results are shown in Table I. There 
is an obvious uncertainty in the given numbers of Table 
I, that are related to our approximate theoretical fitting 
and experimental data fluctuations. However, based on 
such results, the behavior of n(t r ) can be approximately 
described by a linear or exponential decreasing function. 
Here we consider the following simple functional time- rise 
dependences of k: 



K(t r ) = 2300 exp (-0.01 x ut r ) 
K(t r ) = 1900 - 7ujt r . 



or 



(4) 
(5) 



With such decreasing functional time-rise dependence 
of k, we have verified, as shown in Fig. [3J that we have 
the qualitative behavior observed in the experimental re- 
sults given in Ref. |14| of the remaining number of atoms 
N r versus t r . We have considered several values of the 
hold time th- For very small hold time, we can also verify 
the same experimental results that presents a minimum 
of N r as a function of t r . With a linear functional t r de- 
pendence of k, the same qualitative behavior can also be 
reproduced; which differs quantitatively from the expo- 
nential behavior that we show. So, we conclude that the 
time dependent mean field approach can describe all the 
experimental data, if the three-body recombination coef- 
ficient depends on the rise time in this short time scale. 
Such dependence on the rise time can be explored, once 
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FIG. 3: Remaining number of atoms in 85 Rb BEC versus the 
scaled rise time (factor 1/4) of the applied magnetic pulse, for 
some hold times, with hold field 156.7 G (2300ao) and initial 
number Ao = 16500. In the upper frame, we consider k given 
by the exponential dependence of Eq. @; in the lower frame, 
we consider k given by Eq. 0. 

the uncertainties in the experimental results are reduced 
and by considering a better fitting of data (improving 
the values given in Table I). In our interpretation, the 
higher values of k for smaller values of t r are indicating 
the coherent formation of another species (dimers) in the 



condensate. In summary, we report in this work indica- 
tions based on our calculations that a recent experiment 
realized in JILA [14J is evidencing the coherent formation 
of dimers from inelastic collisions in 85 Rb Bose-Einstein 
condensates (BEC). We have solved numerically the non- 
conservative Gross-Pitaevskii equation in spherical sym- 
metry, for condensed systems with very large repulsive 
two-body interaction, varying in time, due to application 
of magnetic field pulses, according to Eq. QJ. According 
to the theoretical predictions of three-body recombina- 
tion rates, we used a dissipation parameter proportional 
to the quartic power of the scattering length a(t) and 
the observed experimental pulse shapes to calculate the 
time evolution of the remaining number of atoms in the 
condensate. We studied this observable as a function of 
the hold and rise times. The experimental results of N r 
versus th can be described in the mean-field approach 
only with very large K, when the rise time is small. This 
indicates a coherent formation of dimers in BEC. For 
longer pulses, when the coherent dimer formation tends 
to disappear, we found that k approaches the maximum 
value given by the theoretical predictions for large scat- 
tering lengths, k ~ 70 |22]- We parameterize this 
surprising behavior of the three-body recombination rate 
considering n = K[t r ) and so it was possible to describe 
the property of a lower dissipation for longer pulses. Fi- 
nally, it is natural to see a burst of relatively hot atoms 
and dimers (carrying the excess of dimer binding energy) 
accompanying a significant loss from the condensate for 
short pulses when the coherent dimer formation occurs. 
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